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We study the stabiUty region of the topological superfluid phase in a trapped two-dimensional 
polarized Fermi gas with spin-orbit coupling and across a BCS-BEC crossover. Due to the competi- 
tion between polarization, pairing interaction and spin-orbit coupling, the Fermi gas typically phase 
separates in the trap. Employing a mean field approach that guarantees the ground state solution, 
we systematically study the structure of the phase separation and investigate in detail the optimal 
parameter region for the preparation of the topologically non-trivial superfluid phase. We then cal- 
culate the momentum space density distribution of the topological superfluid state and demonstrate 
that the existence of the phase leaves a unique signature in the trap integrated momentum space 
density distribution which can survive the time-of-flight imaging process. 
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I. INTRODUCTION 

The study of non-Abelian topological order has at- 
tracted a great amount of interest recently, due to the 
potential applications in fault-tolerant quantum compu- 
tation [2 12]. In addition to systems with intrinsic chi- 
ral p-wave pairing order, e.g. fractional quantum hall 
systems [3^5], chiral p-wave superconductors O [7], p- 
wave superfluidity in ultracold fermions [SHllj etc., it 
has been shown that a topologically non-trivial super- 
fluid phase that supports non-Abelian excitations can be 
induced from an underlying s-wave superfluidity. An ex- 
ample is the semiconductor/superconductor heterostruc- 
tures with spin-orbit coupling (SOC), s-wave pairing su- 
perfluidity and an external Zeeman field [T^lfTS] . With 
the rapidly developing toolbox available for the quan- 
tum control of ultracold atomic systems, the elements 
above can now be implemented experimentally in ultra- 
cold Fermi gases. Importantly, the spin-orbit coupling in 
ultracold Bose gases has been made possible by the recent 
experimental achievement of synthetic gauge field in ul- 
tracold atoms and has generated considerable amount of 
theoretical interests [171 - 120] . On the other hand, the pair- 
ing superfluidity in an ultracold Fermi gas and the quan- 
tum phase transition in a polarized Fermi gas have been 
investigated extensively during the past decade [3TH1H]- 
With clean environment and highly tunable parameters, 
ultracold Fermi gases may serve as an ideal platform for 
the observation of topological superfluidity and for the 
study of the interesting physics therein. In particular, it 
has been suggested that the Majorana zero modes, which 
have eluded experimental observation for decades, may 
be detected in an ultracold atomic system with s-wave 
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interactions [291 - I32] . 

Spin-orbit coupled Fermi gas has been under intensive 
theoretical study recently [311 1551 - H7] . For an unpolar- 
ized Fermi gas near a wide Feshbach resonance, the SOC 
has been found to result in a BCS-BEC type crossover 
even on the BCS side of the resonance [35l [37] . Further- 
more, it has been suggested that the topological super- 
fluid (TSF) phase can be stabilized in a spin polarized 
Fermi gas in the presence of SOC jSTHMl |3i JMl SI]- 
For a polarized Fermi gas without SOC, phase separation 
takes place near a wide Feshbach resonance in a uniform 
gas, due to the competition between population imbal- 
ance and pairing interactions [281 148H50] . With the in- 
troduction of SOC, the phase separation develops a rich 
structure involving the topologically non-trivial super- 
fluid state m] . For a polarized Fermi gas in an ex- 
ternal trap, which is always the case in experiments, the 
various phases naturally separate in real space, with dif- 
ferent phases occurring at different places in the trapping 
potential [HIl [12 [H] . The important questions here are 
whether the topological superfluid is stable in a trapping 
potential and what the detailed structure of the phases 
in a trap is. The phase structure involving the topo- 
logically non-trivial superfluid phase in a trapped three- 
dimensional (3D) polarized Fermi gas with SOC has been 
examined previously, where it has been found that two 
distinct types of topological superfluid may be stabilized 
[39l 141] . In this work, we focus on the phase structure 
of a spin-orbit coupled two-dimensional (2D) polarized 
Fermi gas near a wide Feshbach resonance in a trapping 
potential. 

We study the system using a BCS-type mean field 
theory at zero temperature. Due to the existence of 
metastable or unstable solutions of the gap equation that 
are typical in the presence of population imbalance, we 
directly minimize the thermodynamic potential [28l 141] . 
In contrast to the 3D case where there are two distinct 
topologically non-trivial phases with either two or four 
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gapless points in the quasi-particle excitation spectrum, 
in 2D we find that there is only one topologically non- 
trivial superfluid phase which is always protected by an 
excitation gap away from its phase boundary against the 
conventional superfluid state (SF). This agrees with the 
previous calculations [T2j [16] . As has been shown in Rcf . 
[12j . when a vortex is created in this phase, a Majorana 
zero mode can be found at the center of it. Due to the 
competition between population imbalance and pairing, 
the various phases appear at different locations in an ex- 
ternal trapping potential. We then map out the phase 
diagrams illustrating the structure of the phase separa- 
tion in typical trapping potentials across a wide Feshbach 
resonance under the local density approximation (LDA). 
From our zero temperature mean field calculations, we 
investigate the parameter regions for the existence of a 
stable topological superfluid phase (TSF) in the trap. 
To characterize the properties of different phases, we cal- 
culate the momentum space density distribution for the 
various phases. Notably, the momentum space density 
distribution for the minority spin features a dip near the 
origin in the topological superfluid state, which is unique 
among the phases that we consider. We further demon- 
strate that for appropriate parameters such that the cen- 
ter of the trap is occupied by the TSF state, this signa- 
ture dip of the momentum distribution can survive the 
process of timc-of-flight imaging and thus may serve as 
an unambiguous signal for the TSF. Note that the aim of 
the current work is to provide a qualitatively correct gen- 
eral picture of the phase structure in a trap, as mean field 
theory is not quantitatively reliable near a wide Feshbach 
resonance in 2D due to large fluctuations. We have also 
not considered the possibility of pairing states with non- 
zero center of mass momentum in our mean field theory. 

The paper is organized as follows: in Sec. II, we write 
down the model Hamiltonian and outline the mean field 
theory that we adopt; in Sec. Ill, we analyze the minima 
of the thermodynamic potential and discuss the compe- 
tition between the various phases; we then map out the 
phase diagram of a homogeneous system using the chem- 
ical potential and SOC strength as parameters in Sec. 
IV, which provides us with valuable information regard- 
ing the phase separation in a potential trap; in Sec. V, 
we derive a set of universal dimensionless number equa- 
tions, with which we study the real space distribution of 
particle density and pairing gap in a trap; we investigate 
the momentum distribution of the phases in Sec. VI, and 
identify the signature of the TSF state in the momentum 
distribution; finally, we summarize in Sec. VII. 



II. FORMALISM 

We first consider the model Hamiltonian for a uniform 
2D polarized Fermi gas with Rashba spin-orbit coupling 
near a wide Feshbach resonance. The Hamiltonian can 
be expressed as a sum of three parts: the unperturbed 
Hamiltonian Hq, the SOC Hamiltonian iJgoc and the in- 



teraction Hamiltonian Hint [221 [Ml dS] '■ 

H — ^ ^o-A^a = Ho + Hsoc + Hint 
a 

k,(T k 

+ V ^ «Lt"-k,;"-k',;"k',t, (1) 
k,k' 

where the kinetic energy ek = f?k'^ j (2m), /i^ is the chem- 
ical potential for atoms with spin a = {t, i}, denotes 
the total number of particles with spin u, a\^jj(a\^ ^) an- 
nihilates (creates) a fermion with momentum k and spin 
CT, V is the quantization area in 2D, and H.C. stands for 
Hermitian conjugate. The Rashba spin-orbit coupling 
strength a can be tuned via parameters of the gauge- 
field generating lasers [18], while (/^k = arg (fc^; -I- ifcy). 
In writing the interaction Hamiltonian -ffjnt, we assume 
an s-wave contact interaction between the two fermion 
species, with the bare interaction rate \J renormalized 
following the standard relation in two dimensions [51] : 

k 

Here, Ei, > is the binding energy of the two-body bound 
state in two dimensions without SOC. By tuning through 
a Feshbach resonance from a high- field BCS side, Ei, in- 
creases from zero and becomes large in the low-field BEC 
limit. Therefore, we use the variation of E^ to represent 
the BCS-BEC crossover in the following discussions. One 
should notice that the Ei, we use in this manuscript is not 
the binding energy of the two-body bound state in the 
presence of SOC. 

The non-interacting Hamiltonian Hq + i?soc can be di- 
agonalized in the helicity basis: 

ak,t = ~^ ' 

ak,i (ak,+ - flk,-) , (4) 

where ak,± (ajk-i-) are the annihilation (creation) oper- 
ators for the dressed spin states with different helicities 
(±). Under this basis, the interaction Hamiltonian can 
be written as 

k,k' 

X e""^''' (a_k',+ak\+ - a-k'.-ftk' -) • (5) 
Taking the pairing mean field 

A = ^ ^ (e^"^'^ (a_k,-i-ak,+ - a_k,-ak -)) 

k 

= UY^ («-k,;ak,t) > (6) 

k 
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the mean field Hamiltonian becomes 

k,A=± 

"A 



E 



k,A=± 

"^■^ (a-k,+ak,+ - a-k -ak.-) + h.c. 
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V- 
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where we have defined the chemical potentials — 
(^^ + and h = ^-f - and ^± = ^k ± ak 

with = £k ~ The mean field Hamiltonian is 
quadratic and can be diagonalized in the helicity basis: 

|ak,+,alk,+ '«k,-,aLk-} : 



k,A=± 



k,A=± 



U 



(8) 



Here, Ofk.o- (o^k o-) annihilation (creation) opera- 

tor for the quasi-particles. The quasi-particle excitation 
spectra take the form 



k,± 



i1 + a^k? 



|Ap + ^±2£;o, (9) 



where = y/ih"^/^ + a^k^)^l + h^\A\^/4. From this 
dispersion relation, we see that for finite pairing gap and 
SOC strengths, the only possible gapless point in 2D lies 
in the -Ek,- branch at A: = 0. This takes place when 
h/2 — ^/JJ?~+~A^ . As the chemical potential imbalance 
h increases across this point, the excitation gap first van- 
ishes and then opens up again, and the system enters a 
topologically non-trivial superfluid phase [HI . This 
is in contrast to the 3D case, where two distinct topolog- 
ically non-trivial phases exist, with two or four gapless 
points in the quasi-particle excitation spectrum 39, 41]. 

In this manuscript, we consider only the zero temper- 
ature case and write down the thermodynamic potential 
as 
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where j3 ~ l/ksT and ks is the Boltzmann constant. 
Considering the cxtrema condition of the thermody- 
namic potential dfl/dA — and the number constraints 
Uc = (— l/V)9^2/9/icr, we get the gap and the number 
equations, respectively: 
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where 5^ = —S^ = — 1, and we have taken A to be real 
for simplicity. The ground state of the system at zero 
temperature is given by the global minimum of the ther- 
modynamic potential in Eq. ( 10 1 under the number con- 
straints Eq. (12). For a uniform gas, one has to con- integration: N„ = J (Prna{v), where the local density 



center ^ and the chemical potential imbalance h are re- 
lated to the total particle number N ~ + and the 
polarization P = {N^ — N^) /N. The total particle num- 
ber for each spin species can be determined from a trap 



sider explicitly the possibility of phase separation and 
introduce a mixing coefficient in order to get the correct 
ground state. In an external trapping potential, the var- 
ious phases naturally separate in real space SU [SD] . 

Next, we focus on the phase separation in the pres- 
ence of an external trapping potential V(r), due to its 
experimental relevance. Assuming the potential to be 
slowly varying and taking the local density approxima- 
tion (LDA) , we can write the chemical potentials at each 
spatial location r as: yU^(r) = fij. + h/2, /Lt^(r) = yUp — /i/2, 
and /ir = — V{r), where the chemical potential at trap 



no-(r) can be calculated from Eq. (12 1 with /i replaced 



by /ir, and with the local pairing order parameter A(r) 
determined from the global minimum of the thermody- 
namic potential at each spatial location r. Without loss 
of generality, we assume > throughout this work 
such that h and P are both positive. 
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FIG. 1: (Color online) Typical shapes of thermodynamic po- 
tential as a function of the pairing order parameter A with 
various SOC strength akh/h (left column) and chemical po- 
tential ^/h (right column). The two-body binding energy is 
chosen as Eb/h = 0.5. The parameters of the subplots are: 
(a)/i//i = 0.2, akh/h = 0.1, ih)n/h = 0.2, akh/h = 0.35, 
(c)/i//i = 0.2, akh/h = 0.6, {d)ii/h = 0.2, akh/h = 0.8, 
{e)fi/h = 0.4, akh/h = 0.3, (f)p/ft = 0.24,afch//i = 0.3, 
(g)/i//i = 0.1, akh/h = 0.3, (h)/i//i = -0.2, akh/h = 0.3. The 
chemical potential difference h is taken to be the energy unit, 
while the unit of momentum kh is defined as ft^fe^/(2m) = h. 



III. THERMODYNAMIC POTENTIAL IN THE 
PRESENCE OF SOC 

In Fig. [l] we show the behavior of the thermodynamic 
potential for a set of typical parameters. Similar to the 
case of a polarized Fermi gas without SOC, the competi- 
tion between polarization and pairing leads to a double- 
well structure in the thermodynamic potential in certain 
parameter regions [c.f. Fig. TJa)(b)(e)(f)(g)]. As a conse- 
quence, two distinct gapped phases can be present in the 
phase diagram, which are separated by a quantum phase 
transition. Specifically, if one of the two gapped states is 
a conventional superfluid with h/2 < \//z^~-j-~A^, while 
the other a topological superfluid with h/2 > \/ fJ-^ + A^, 
there must be a first order phase transition between SF 
and TSF phases as the parameters are tuned so that the 
ground state of the system changes from one local min- 
imum to the other. If both of the pairing orders are 
conventional superfluid with the same symmetry, there 
can only be a first-order-like crossover [16j. Due to the 
non-monotonic behavior of the thermodynamic potential, 
the solution of the gap equation may also correspond to 
metastable or unstable states, in addition to the ground 
state. To avoid this complication, we directly minimize 
the thermodynamic potential to ensure that the ground 
state is reached. 
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FIG. 2: (Color online) Dependence of the pairing gap on SOC 
strength, (a) Evolution of the gap equation with increas- 
ing SOC strengths. Form the upmost curve to the lowest: 
akh/h = 0.8,0.5,0.3,0.2, respectively, (inset) Enlarged view. 
G is related to the derivative of the thermodynamic poten- 
tial with respect to A, as defined in the text, (b) Scaling of 
the SOC-induced pairing order with SOC strength, (inset) 
Semi-logarithm plot of the scaling relation. For both plots, 
Eb/h = 0.5, ii/h = 0.1. 



In the absence of SOC, when the polarization becomes 
large enough, a phase transition occurs and brings the 
system from a superfluid phase to a normal phase |28) . 
However, an arbitrarily small SOC will change this pic- 
ture and introduces novel type of phases and phase tran- 
sitions to the system. In fact, when A = 0, a singularity 
exists in the integrand of the gap equation ( 11 ) over con- 



siderably large parameter regions. So long as this singu- 
larity exists, the gap equation always has at least one fl- 
nite solution regardless of the SOC strength a and chemi- 
cal potential combinations (//, h). In order to understand 
this picture, we show in Fig. [2]ja) the behavior of the 
function G = (— l/2A)(?ri/(9A, which is proportional to 
the left-hand side of the gap equation (11). In the pres- 



ence of the singularity, for arbitrary SOC strength, the 
function G is always diverging as A — > 0, and tends to 
large negative values as A — )■ oo. Therefore, there is at 
least one solution to the gap equation under these con- 
ditions, indicating the presence of gapped phases, as has 
been pointed out previously in Ref |16j . Further analysis 
shows that one of the gapped phases is the global mini- 
mum of the thermodynamic potential. This observation 
shows that superfluidity can survive arbitrary polariza- 
tion, provided that an SOC is introduced. In Fig. [2]^b), 
we present the pairing gap A of the ground state as a 
function of SOC strength a. The numerical result sug- 
gests that the pairing gap decreases super-exponentially 
as a approaches zero. 

The singularity responsible for the divergence of G 
comes from the terms proportional to E^\, which di- 
verges at A = 0. As Sk,- = ll^kl - ^/^^P^WTh?Jl\ 
at A = 0, the solutions of the equation |^k| = 

a^k/^ + /i^/4 give the singularity points in momen- 
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FIG. 3: (Color online) The phase diagram in the a-^ plane 
with the binding energy Eb/h = 0.5. The first order phase 
transition is shown in red solid curve while the second order 
phase transition in dash-dotted black curve. The thin dashed 
curve in the TSF region marks the A/h — 10~^ threshold. 
The chemical potential difference h is taken to be the en- 
ergy unit, while the unit of momentum kh is defined through 
h'kl/{2m) = h. 
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FIG. 4: (Color online) The phase diagram in the a-^L plane 
with various binding energies E^. {sL)Et/h = 0.2, {h)Et/h = 
0.65, {c)Et/h = 0.95, {d)Eb/h = 1.2. The thin dashed curves 
in (a-c) mark the A/h — 10"'' threshold. The chemical po- 
tential h is taken to be the energy unit, while the unit of 
momentum kh is defined through fi^fc^/(2m) — h. 



turn space. It is easy to see that the equation above 
does not have real-valued solutions under the conditions 
< -{a^ /i2)/(4a2) or < min(-/i/2, -aV2). This 
suggests that the ground states corresponding to Fig. 
[lja)(g)(h) are superfluid phases with small but finite 
pairing gap. 



IV. PHASE DIAGRAM IN THE a-^i PLANE 

From the discussions in the previous section, we see 
that the presence of SOC can lead to a rich structure of 
phases in a trapping potential. As a first step to under- 
stand the spatial distribution of the various phases in the 
trap, we consider in this section a homogeneous system 
and investigate the phase diagram as a function of (a, fi) 
for given Ei, and h. Under LDA while assuming both 
spin species experience the same harmonic potential, a 
downward vertical line in such a phase diagram repre- 
sents a trajectory from a trap center to its edge, with the 
chemical potential at the trap center fixed by that at the 
starting point of the line. To this end, we only need to 
minimize the thermodynamic potential in Eq. ( 10 1 for 



given SOC strength a and chemical potential difference 
h while sweeping the chemical potential /x. 

In Fig. [3) we show a typical phase diagram in the 
a-/i plane for Ef,/h = 0.5. Notice that there is only 
one topologically non-trivial superfluid phase, which is 
clearly different from the 3D case as discussed before. 
The TSF phase is separated from the conventional super- 
fluid phase by two kinds of phase boundaries. The solid 
curve in Fig. [3] represents a first order phase boundary, 
along which the states corresponding to the two local 
minima of the double well structure in the thermody- 
namic potential are degenerate in energy. Compared to 
the 3D case, the first order phase boundary is dramati- 
cally extended. The other kind of TSF-SF phase bound- 
ary is of second order, given by h/2 = v'/i^ A2, along 
which the pairing gap remains finite and the excitation 
gap vanishes. 

To determine the phase boundary for A = 0, we 
need to examine the existence of divergence in the gap 
equation as A approaches zero. As discussed in the 
previous section, the singularities go away when /i < 
-(a^ + h?)/{'ia^) or ^ < m\n{~h/2,~a^ /2). The phase 
boundary of the superfluid phases with A = can be 
calculated by solving the gap equation in these param- 
eter regions. We note that the maximum value of the 
chemical potential satisfying these relations is —h/2, be- 
low which the chemical potential of both spin species ^„ 
become negative. Hence there will not be a phase bound- 
ary between a superfluid phase (SF or TSF) and a normal 
phase. Instead, only phase boundaries between a super- 
fluid state and vacuum (VAC) exist. For the calculations 
above, we always check the thermodynamic potential to 
ensure that states along the phase boundary with A = 
represent ground state solutions. 

According to the phase diagram in Fig. [3| the stabil- 
ity region for the TSF phase appears to be significant. 
Yet this can be misleading for experimental detection. 
In fact, the size of the pairing gap in the TSF phase with 
small SOC strength is typically vanishingly small. This 
can be seen from the dashed curve traversing the TSF 
phase in Fig. |3j which is solved from the gap equation 
by setting A/h = 10"^. To the left of the curve, the pair- 
ing gap A/h < 10~^ and decreases exponentially fast as 
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FIG. 5: (Color online) The distributions of number densities na(f) (a-e) and the order parameter A(f) (f-j) are shown versus 
dimensionless distance from the trap center f = r/R. The parameters for each column are: (a) Eb/Ep ~ 0.32, akp/Ep = 
0.3, h/Ep = 1,P = 0.624; (b) E^/Ef = 0.5, akp/Ep = 0.4, h/Ep = 1, P = 0.287; (c) E^/Ep = 0.5, akp/Ep = 0.7, h/Ep = 1, 
P = 0.188; (d) Eb/Ep = 0.5, akp/Ep = 0.8, h/Ep = 1, P = 0.176; (e) Et/Ep = 0.5, akp/Ep = 0.6, h/Ep = 1.45, P = 0.662. 
The bottom row (k-o) illustrates the shell structure of phase separation. The solid black (dashed red) curves in the density 
subplots represent spin up (down) species. The thin dotted lines in the first two rows illustrate the TSF-SF or the SF-SF 
boundary. The units of energy Ep and of length R are defined in the text, and the unit of density is no = mEp/ {nh'^). 



a approaches zero. The order parameter A only becomes 
significant when a is further increased toward the phase 
boundary between TSF and SF. Given the fluctuations 
in 2D systems at finite temperatures, experimental ob- 
servation of the TSF phase is only possible to the right 
of the dashed curve and with reasonably large pairing 
gap A. 

Figure [3] also provides information regarding the struc- 
ture of the phase separation in a trapping potential. 
When the SOC is small, the Fermi gas will phase sep- 
arate into two regions, a conventional superfluid core 
surrounded by a TSF phase with large spin polarization 
and vanishingly small pairing order. The phase bound- 
ary between them is of first order. As the SOC increases, 
the local minima in the thermodynamic potential corre- 
sponding to the TSF and the SF states move closer as the 
pairing gap of TSF state increases. The two local minima 
merge at a critical end point beyond which the double- 
well structure in the thermodynamic potential disappears 
and the phase boundary between TSF and SF becomes 
second order. Further increasing the SOC, there may be 
a parameter window where the TSF phase appears as a 
ring structure in the trap. Finally, when the SOC is large 
enough, phase separation no longer occurs and the trap 
is filled with a superfluid of rashbons. 

We have also calculated the a-^ phase diagram for a 



homogeneous system with different bound state energies 
Eb (see Fig. ^. Toward the BCS side [Fig. ga)], the 
stability region of the TSF phase increases while the first 
order phase boundary between TSF and SF no longer ex- 
ists. For small SOC and large chemical potential, there 
may exist two different SF phases at the trap center, 
separated by a first-order-like boundary. In this case as 
the symmetries of the two SF phases are the same, the 
boundary is merely a crossover. On the phase diagram, 
this first-order-like crossover boundary ends at a critical 
end point where the two potential wells in the double well 
structure of the thermodynamic potential merge. Im- 
mediately below this first order crossover boundary and 
with small SOC, an SF phase with vanishingly small or- 
der parameter and small polarization appears where the 
chemical potentials of both spin species are positive. This 
corresponds to an SOC induced SF phase out of a nor- 
mal phase with two spin species without SOC. Toward 
BEC side [Fig. gb-d)], the stability region of the TSF 
phase becomes smaller and eventually disappears from 
the phase diagram. The trap is then occupied by super- 
fiuid of rashbons. 
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FIG. 6: (Color online) The phase structure appearing in a 
harmonic trapping potential with the parameters Eb/Ep = 
0.5 and akp/Ep — 0.75. Here, the total polarization P is a 
trapped integrated result as calculated from Eq. (141, and 
f 



r/R is the dimensionless distance from the trap center. 
Notice that the trap is filled with TSF phase as P is above a 
critical value. All phase boundaries shown in this plot are of 
second order. First order phase boundaries show up at smaller 
SOC strengths and/or smaller Et- The units of energy Ep 
and of length R are defined in the text. 
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FIG. 7: (Color online) The density distribution in momen- 
tum space, (a) TSF with small gap: Eb/h = 0.5, fj,/h = 0, 
akf/h = 0.1; (b) TSF with larger gap: Eb/h = 0.5, n/h = 0, 
akf/h = 0.45; (c) SF: Eb/h = 0.5, fi/h = 0, akf/h = 0.8, (d) 
SF with smaU gap: Eb/h = 1.2, fi/h = 0, akf /h = 0.1, (e) 
SF: Eb/h = 1.2, fi/h = 0, akf/h = 0.45, (f) SF: Eb/h = 1.2, 
fj,/h — 0, akf/h — 0.8. The chemical potential h is taken 
to be the energy unit, the unit of momentum kh is defined 
through h^k^/{2m) = h, and the unit of density is defined as 
n„ = kl/{27T). 



Next, we adopt the LDA and explicitly include the 
trapping potential in our calculation. To make our cal- 
culation universal and applicable to systems with any 
total particle number, we derive a dimensionless form 
following Ref. [48]. We take the unit of energy to be 
the Fermi energy (Ep) at the center of a 2D axially sym- 
metric trap for N non-interacting fermions with equal 
population for the two spin species, with Ep — huj^/N 
and u is the trapping frequency. The harmonic trapping 
potential in the dimensionless form can be expressed as 
V{r)/EF = r^/R'^ = P, where R = ^jEp/m is the 
Thomas-Fermi radius in two dimensions |52| . The num- 
ber equation in the dimensionless form then becomes 



1 



d2f[n^(f) -f n;(f)], 



(13) 
(14) 



with dimensionless number density h„ = n^/ni^. Here, 
ria is the number density given by the number equation 
Eq. (12) at position r. The Fermi momentum kp is 

It is 



fc|./(27r). 



defined as Ep — h kp/{2'm), and no 
obvious that the properties of the system only depend 
on the dimensionless parameters {Eb/ Ep, akp/ Ep, P}. 

Solving the dimensionless equations above, we get the 
typical phase structure in a trapping potential with a var- 
ious sets of parameters. Note that for simplicity, we first 
choose an appropriate chemical potential difference, e.g. 



h = Ep, and solve for the chemical potential fi at the 
center of the trap from Eq. ( 13 ) with fixed SOC strength 
a and Ei,. The polarization P can then be calculated 
from Eq. ( [T4| ). The resulting shell structures are shown 
in Fig. [5] The topological superfluid phase typically ap- 
pears toward the edge of the trap or as a ring between 
two SF phases, in agreement with our previous discus- 
sions based on the phase diagram of Fig. |3] Notably, 
there are parameter regions where the TSF state can oc- 
cupy the entire trap. This is shown in the right-most 
column in Fig. [Sj and corresponds to a vertical line in 
the a-fi plane phase diagram with the starting point in 
the TSF phase. 

To better understand the phase structures in a trap, we 
plot in Fig. [6] the zero temperature phase boundaries in 
a harmonic trapping potential in the P-f plane, where P 
is the trap-integrated total polarization calculated from 
Eq. (14), and r = r/Ris the dimensionless distance from 



the trap center. Notice that the TSF phase occupies the 
entire trap only when the total polarization exceeds a 
critical value (~ 0.69 for the chosen parameters in Fig. 
[g]). As we will show in the next section, this regime pro- 
vides an ideal setup for the detection of the TSF state 
in a trapped gas. When P decreases from this value, the 
conventional SF phase will emerge from the trap center, 
gradually extend to the trap edge, and eventually occupy 
the entire trap for small polarization case. 
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VI. MOMENTUM DISTRIBUTION AND THE 
SIGNATURE OF THE TOPOLOGICAL 
SUPERFLUID PHASE 

To characterize the properties of the various phases 
in the phase diagram, we calculate their respective mo- 
mentum distribution (see Fig. [7|), which is given by the 
summand in the number equations (12). In Fig. [t] we 



(a) 



show the momentum distribution of a homogeneous sys- 
tem with various parameters. In the first row of Fig. [Tj 
the binding energy is set as Ei,/h = 0.5 with increasing 
SOC strength. In the second row, a similar evolution 
with a is shown but with a binding energy Eiy/h = 1.2, 
more toward the BEC regime. It is apparent that the 
momentum distribution in a TSF phase [c.f. Fig. [TJa-b)] 
is drastically different from that in an SF phase. In par- 
ticular, the momentum distribution of the minority spin 
in the TSF phase features a dip near zero momentum, 
which can be explained by the observation that rik,^ = 
at fc = 0, where lis the summand in the corresponding 
number equation ( |12[ ). 

As this dip in the momentum distribution is unique 
to the TSF phase, one may think of using it as a sig- 
nature for the experimental detection of the TSF phase. 
To measure the momentum distribution experimentally, 
a commonly used practice is the time-of-flight imaging 
technique, which involves a ballistic expansion of the gas 
after suddenly switching off the trapping potential. As 
there are typically several different phases in the trap- 
ping potential, the observed momentum distribution is 
usually a trap-integrated distribution which includes the 
contribution from all the phases in the trap. In this case, 
the signal of the topological superfluid is washed out and 
cannot be detected. 

One possible way to overcome this difficulty is to pre- 
pare the system in an appropriate parameter region such 
that the center of the trap is filled with the TSF state. 
An example of this is demonstrated in Fig. [sj^b). The 
corresponding trap-integrated momentum distribution is 
shown in Fig. [sjja), where the signature of the TSF state 
apparently survives the trap integration. In comparison, 
we show in Fig. |8][c)(d) similar calculations for the case 
with an SF core surrounded by the TSF phase. As is 
clear from Fig. [sjjc), the signature dip for the TSF state 
can no longer be observed in the trap-integrated momen- 
tum distribution. This suggests that the existence of the 
dip can serve as a signature for the existence of the TSF 
phase if the momentum distribution of the minority spin 
species can be detected. 
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FIG. 8: (Color online) (left column)The trap-integrated den- 
sity distribution in momentum space Uk- (right column) The 
number density distribution for both spin up (blaclc solid 
curves) and spin down (dashed red curves) atoms in a trap- 
ping potential. The insets show the distribution of the order 
parameter A(f). The dashed line in (d) illustrates the TSF- 
SF phase boundary. The parameters are: {a.){b)Eb/ Ef ~ 0.2, 
akp/Ep = 1.25, P = 0.82; (c)(d) Et/Ep = 0.2, akp/Ep = 1, 
P = 0.65. The energy unit Ep is defined in the text and the 
unit in momentum space /cf is related to the unit of energy 
as h'^kp/{2m) = Ep- The unit of density is defined through 
no = fc|/(27r). 



rameters. Compared to the 3D case, we find dramatically 
increased first order phase boundary between the SF core 
and the TSF phase that surrounds it, which makes it ob- 
servable in experiments from the density distributions of 
the spin species. We then develop a universal scheme 
for the characterization of a trapped gas. The resulting 
phase and density distributions are therefore independent 
of the trapping geometry and the total particle number, 
and are determined by a set of dimensionless parameters. 
We explicitly calculate the density and momentum dis- 
tribution of the gas in a trapping potential. Importantly, 
we find a parameter region where the trap is occupied 
by the TSF phase only. In this regime, the characteristic 
signature of the TSF state in the momentum distribution 
can survive the trap integration, rendering the signal de- 
tectable in a time-of-flight imaging process. 
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